Collaborators: Tomo Tanaka, Hari Doodhi
Notebooks: 1
In this document we present a simple computer model for error correction in kinetochore-microtubule interaction during bi-orientation. The configuration is very simple, perhaps simplistic, but it can test the underlying mechanisms. We focus on two situations:
In our model there are two kinetochores, left and right. Each of them can be in one of the following attachment states:
Dual attachments (ambilink, synlink, biend) are present in the second model, when the lateral attachment is formed on top of the existing end-on attachment. In addition to the connection state, we also record which spindle pole (left or right) the kinetochore is attached to. Therefore, the state of the system can be briefly denoted as, e.g.:
lateral-R endon-R
The left KT has a lateral attachment to a microtubule extending from the right spindle. The right kinetochore has an end-on attachment to a microtubule extending from the right spindle.
Transition between states are governed by a Poisson law. That is, at any given time the probability of transition from, e.g., lateral to end-on attachment is constant. Each transition happens randomly a rate, \(R\), expressed in events per minute.
The model allows for other probability prescriptions, for example, the transition time can be expressed by a Gaussian probability.
The simulation is a sequence of events. The time of the next event is randomly drawn based on the current event and the probability law assumed. The events are generated for both kinetochores concurrently, until the final stable configuration is established (endon-L endon-R).
Here we assume that the end-on detachment and formation of the lateral attachment are independent. For illustration purposes, we show how the model works on the left kinetochore, but obviously, it works on both. Also, microtubules can reach kinetochores at the opposite side. The sequence of events in this model is as follows.
Let’s start with an unattached microtubule.
A lateral attachment is created at formation rate \(R_\mathrm{form}\).
Next, the lateral attachment is converted into the end-on attachment, at conversion rate \(R_\mathrm{conv}\):
Finally, the microtubule is detached at detachment rate \(R_\mathrm{det}\):
In this model we assume that the end-on detachment happens only upon formation of the lateral attachment. Hence, we never observe a detached state. This is the sequence of events. Now, we start with an end-on attachment:
Then, a MT from the opposite spindle can make a lateral takeover, at replacement rate \(R_\mathrm{rep}\). This is the same process as “formation” in model 1, but we call it “replacement” here, for distinction. For a moment we have an “ambilink” attachment.
The ambilink attachment is unstable and leads to the lateral attachment by knocking off the end-on MT at a rate \(R_\mathrm{ko}\):
As before, the lateral attachment is converted into end-on at conversion rate \(R_\mathrm{conv}\).
The lateral takeover switches the side of the MT. E.g., it starts with an endon-L configuration and ends up with endon-R configuration.
When the lateral takeover is attempted from the same side, observations suggest that it fails to knock-off the existing end-on attachment. Instead, a “synlink” attachment is formed (at replacement rate, \(R_\mathrm{rep}\)):
This works like a normal lateral attachment and is converted into endon at conversion rate \(R_\mathrm{conv}\), just like above, except that at the end of the process a “biend” (bi-end-on) attachment is formed.
This configuration is unstable and the new MT (the one that was trying to take over) is the detached at the knock-off rate \(R_\mathrm{ko}\):
The failed lateral takeover does not change the side of the MT.
Note: both opposite-side and same-side lateral takeover take the same amount of time. This is because the opposite-side takeover is a sequence of events with rates \(R_\mathrm{rep}\), \(R_\mathrm{ko}\) and \(R_\mathrm{conv}\), while the same-side takeover is a sequence of events with rates \(R_\mathrm{rep}\), \(R_\mathrm{conv}\) and \(R_\mathrm{ko}\). The only difference is the switched order of the last two events, but the total time (on average) will be the same.
We estimate conversion rate (from lateral to end-on) from our old simulations of kinetochore-microtubule interaction. We use the latest results from the default parameter model in this file:
## [1] "/home/mgierlinski/projects/kinetube/varcap_grid/varcap8_slr2_short0.2_long0.5_ktd0.1_mtd0.03_rf0.3_resp0.6_inimt10_tau0.5_beta0.7/kinetube_all_kt.dat"
It contains several statistics for each kinetochore, over many simulations. Among them is the total sliding time. This is not exactly what we want, as this is the accumulated time, so includes KTs that were detached and re-attached. But the simulation also recorded the total drop-off time, that is the time the KT spend while being knocked-off before being re-attached. We select only data where total drop-off time is zero, which corresponds to KTs that have not been knocked-off. The mean sliding time from this file is 0.67 minutes, which gives the conversion rate of 1.5 per min.
I modified the old simulation script to record conversion times, from onset of sliding until end-on conversion. This excludes pole capture. The new script is
/cluster/gjb_lab/mgierlinski/projects/kinetube/scripts/simconv.pl
and to run it on the cluster we used the following
qsub -q c6100.q,c6145.q -R y -V -cwd -t 1-100 ./scripts/simconv.pl -n -nsim=30 -delta=0.001 -poletransfer -ktexclusion -mtshrinkpole=0 -initialdiffusion=0.1 -maxtime=30 -alpha=1.57 -ktradrate=2 -mtdiffusion=0.03 -diffusion=0.1 -rfrate=0.3 -ktresprob=0.6 -inimt=10 -caprads=0.2 -capradl=0.5 -cumuldisttau=0.5 -beta=0.7 -dir=convrate -root=convrate
Technical note: I had to limit the queues, because all the jobs sent to c6320.q dies with “No such file or directory”. I suspect Perl configuration is different on these nodes.
Here is the distribution of individual conversion times (the same distribution in linear of logarithmic scale):
The mean conversion time is 0.488 min, which gives conversion rate of 2.05 per min.
First, I test my internal function for accuracy. Using all default parameters. The red curves show the theoretical exponential distribution.
Here is a table with mean and SE for these:
| event | mean duration | SE |
|---|---|---|
| formation | 0.4879 | 0.0049 |
| conversion | 0.5035 | 0.0051 |
| detachment | 1.0125 | 0.0102 |
| knockoff | 0.0000 | 0.0000 |
| replacement | 0.5006 | 0.0050 |
We undertake a few simple tests to confirm that the model is working correctly. Specifically, we build the distribution of event duration over many simulation runs. The mean duration should be as in the input parameters. We set the following parameters: \(R_\mathrm{form} = 0.5~\mathrm{min}^{-1}\), \(R_\mathrm{conv} = 2~\mathrm{min}^{-1}\) and \(R_\mathrm{det} = 1~\mathrm{min}^{-1}\).
Here are measured distributions of event duration versus theoretical exponential distributions.
And these are mean of these distributions, which should correspond to 1 / rate (that is we expect means of conversion = 0.5, detachment = 1 and formation = 2):
| event | mean duration | SE |
|---|---|---|
| conversion | 0.5030 | 0.0019 |
| detachment | 0.9951 | 0.0033 |
| formation | 1.9941 | 0.0075 |
There is a good agreement between generated and theoretical means.
The same for the “M2” model:
and the mean times:
| event | mean duration | SE |
|---|---|---|
| conversion | 0.5004 | 0.0032 |
| knockoff | 0.0000 | 0.0000 |
| replacement | 1.9961 | 0.0094 |
They are reasonably close to the input parameters. In the second test the “knockoff” is very short: this corresponds to infinitesimal delay of detachment after lateral attachment. When we set a knockoff rate of 1 per minute, we get:
| event | mean duration | SE |
|---|---|---|
| conversion | 0.4980 | 0.0027 |
| knockoff | 1.0103 | 0.0054 |
| replacement | 2.0186 | 0.0086 |
Everything is in good agreement.
We start in an error configuration, with both kinetochores connected to the same spindle:
Here is an example of a simulation, showing the sequence of states for each kinetochore. The “L” or “R” in each state denotes the spindle to which the kinetochore is connected.
| time | A | B |
|---|---|---|
| 0.00 | endon-L | endon-L |
| 0.11 | — | endon-L |
| 0.23 | lateral-L | endon-L |
| 0.27 | lateral-L | — |
| 0.44 | lateral-L | lateral-L |
| 1.44 | lateral-L | endon-L |
| 4.55 | lateral-L | — |
| 4.94 | lateral-L | lateral-L |
| 5.52 | endon-L | lateral-L |
| 5.57 | — | lateral-L |
| 5.92 | lateral-L | lateral-L |
| 6.17 | endon-L | lateral-L |
| 6.33 | — | lateral-L |
| 6.64 | lateral-L | lateral-L |
| 7.01 | lateral-L | endon-L |
| 7.66 | lateral-L | — |
| 7.78 | lateral-L | lateral-R |
| 8.54 | endon-L | lateral-R |
| 8.99 | — | lateral-R |
| 9.05 | lateral-L | lateral-R |
| 10.93 | endon-L | lateral-R |
| 11.19 | endon-L | endon-R |
Here is a graphical representation of this simulation run:
The final configuration is the correct “endon-L endon-R” (“endon-R endon-L” is also considered final).
Another example, this time for model 2. Knockoff rate is high (\(10\)), hence “ambilink” and “biend” states are short.
| time | A | B |
|---|---|---|
| 0.00 | endon-L | endon-L |
| 0.37 | synlink-L | endon-L |
| 0.43 | biend-L | endon-L |
| 0.55 | endon-L | endon-L |
| 0.91 | endon-L | synlink-L |
| 1.04 | endon-L | biend-L |
| 1.11 | endon-L | endon-L |
| 1.47 | endon-L | synlink-L |
| 3.28 | synlink-L | synlink-L |
| 4.91 | biend-L | synlink-L |
| 5.00 | endon-L | synlink-L |
| 5.46 | synlink-L | synlink-L |
| 7.13 | biend-L | synlink-L |
| 7.24 | endon-L | synlink-L |
| 8.06 | endon-L | biend-L |
| 8.08 | endon-L | endon-L |
| 9.23 | endon-L | ambilink-R |
| 9.54 | endon-L | lateral-R |
| 9.72 | synlink-L | lateral-R |
| 10.85 | biend-L | lateral-R |
| 10.88 | endon-L | lateral-R |
| 11.42 | endon-L | endon-R |
Parameter sets:
| parameter | set 1 | set 2 |
|---|---|---|
| formation | 1 | 0.5 |
| conversion | 2 | 2.0 |
| detachment | 1 | 1.0 |
Parameter set 1:
| parameter | set 1 | set 2 |
|---|---|---|
| replacement | 1 | 0.5 |
| knockoff | 10 | 10.0 |
| conversion | 2 | 2.0 |
Here are distributions of simulation time for both models (up and down) and parameter sets (fill colours). We show the same data in linear of logarithmic scale.
Here is the mean and median simulation time for
| model | set | mean | median |
|---|---|---|---|
| M1 | 1 | 6.3 | 4.6 |
| M1 | 2 | 12.3 | 8.8 |
| M2 | 1 | 2.6 | 1.9 |
| M2 | 2 | 3.4 | 2.5 |
We run the simulation across a grid of parameters. This table show the rates used:
| formation rate | 0.125, 0.25, 0.5, 1, 2, 4 |
| conversion rate | 0.125, 0.25, 0.5, 1, 1.5, 2, 4 |
| detachment rate | 0.0625, 0.125, 0.25, 0.5, 1, 2, 4 |
Simply, we start from the default values and vary each parameter by factor two, up and down (except for knockoff). For each set we run 10,000 simulations.
The table below shows results ordered by the mean simulation time. There are 294 combinations of the three parameters.
| formation | conversion | detachment | mean.time | se.time | median.time | median.detached.fraction |
|---|---|---|---|---|---|---|
| 4.00 | 4.00 | 2.00 | 2.02 | 0.02 | 1.50 | 0.01 |
| 4.00 | 4.00 | 4.00 | 2.27 | 0.02 | 1.63 | 0.07 |
| 4.00 | 4.00 | 1.00 | 2.25 | 0.02 | 1.69 | 0.00 |
| 2.00 | 4.00 | 2.00 | 3.12 | 0.03 | 2.27 | 0.10 |
| 2.00 | 4.00 | 1.00 | 3.06 | 0.03 | 2.31 | 0.01 |
| 4.00 | 2.00 | 2.00 | 3.14 | 0.03 | 2.31 | 0.01 |
| 4.00 | 2.00 | 1.00 | 3.08 | 0.03 | 2.32 | 0.00 |
| 4.00 | 4.00 | 0.50 | 3.10 | 0.03 | 2.34 | 0.00 |
| 4.00 | 1.50 | 1.00 | 3.61 | 0.03 | 2.68 | 0.00 |
| 4.00 | 2.00 | 4.00 | 3.96 | 0.04 | 2.82 | 0.04 |
| 2.00 | 4.00 | 0.50 | 3.75 | 0.03 | 2.82 | 0.00 |
| 2.00 | 4.00 | 4.00 | 3.97 | 0.04 | 2.85 | 0.21 |
| 4.00 | 2.00 | 0.50 | 3.81 | 0.03 | 2.88 | 0.00 |
| 4.00 | 1.50 | 2.00 | 4.01 | 0.04 | 2.90 | 0.01 |
| 2.00 | 2.00 | 1.00 | 4.01 | 0.03 | 3.01 | 0.01 |
| 4.00 | 1.50 | 0.50 | 4.27 | 0.04 | 3.20 | 0.00 |
| 2.00 | 2.00 | 2.00 | 4.49 | 0.04 | 3.32 | 0.08 |
| 2.00 | 2.00 | 0.50 | 4.46 | 0.04 | 3.43 | 0.00 |
| 2.00 | 1.50 | 1.00 | 4.69 | 0.04 | 3.51 | 0.01 |
| 4.00 | 1.00 | 1.00 | 5.03 | 0.05 | 3.67 | 0.00 |
| 4.00 | 4.00 | 0.25 | 5.10 | 0.05 | 3.68 | 0.00 |
| 1.00 | 4.00 | 1.00 | 5.11 | 0.05 | 3.71 | 0.12 |
| 2.00 | 1.50 | 0.50 | 4.94 | 0.04 | 3.72 | 0.00 |
| 4.00 | 1.50 | 4.00 | 5.43 | 0.05 | 3.83 | 0.03 |
| 4.00 | 1.00 | 0.50 | 5.17 | 0.04 | 3.85 | 0.00 |
| 1.00 | 4.00 | 0.50 | 5.26 | 0.04 | 3.93 | 0.01 |
| 2.00 | 1.50 | 2.00 | 5.56 | 0.05 | 4.02 | 0.06 |
| 2.00 | 4.00 | 0.25 | 5.66 | 0.05 | 4.15 | 0.00 |
| 4.00 | 2.00 | 0.25 | 5.63 | 0.05 | 4.18 | 0.00 |
| 4.00 | 1.00 | 2.00 | 6.16 | 0.06 | 4.36 | 0.01 |
| 1.00 | 4.00 | 2.00 | 6.11 | 0.06 | 4.38 | 0.27 |
| 1.00 | 2.00 | 0.50 | 6.04 | 0.05 | 4.52 | 0.01 |
| 2.00 | 2.00 | 4.00 | 6.37 | 0.06 | 4.54 | 0.14 |
| 2.00 | 1.00 | 1.00 | 6.21 | 0.05 | 4.56 | 0.01 |
| 4.00 | 1.50 | 0.25 | 6.14 | 0.05 | 4.56 | 0.00 |
| 1.00 | 2.00 | 1.00 | 6.31 | 0.06 | 4.59 | 0.10 |
| 2.00 | 1.00 | 0.50 | 6.06 | 0.05 | 4.62 | 0.00 |
| 2.00 | 2.00 | 0.25 | 6.29 | 0.05 | 4.76 | 0.00 |
| 2.00 | 1.50 | 0.25 | 6.74 | 0.06 | 5.08 | 0.00 |
| 4.00 | 1.00 | 0.25 | 6.83 | 0.06 | 5.14 | 0.00 |
| 1.00 | 1.50 | 0.50 | 6.74 | 0.06 | 5.14 | 0.01 |
| 1.00 | 1.50 | 1.00 | 7.06 | 0.06 | 5.14 | 0.09 |
| 1.00 | 4.00 | 0.25 | 6.98 | 0.06 | 5.18 | 0.00 |
| 1.00 | 2.00 | 2.00 | 7.99 | 0.08 | 5.69 | 0.21 |
| 1.00 | 2.00 | 0.25 | 7.53 | 0.06 | 5.71 | 0.00 |
| 2.00 | 1.50 | 4.00 | 8.00 | 0.08 | 5.75 | 0.11 |
| 2.00 | 1.00 | 2.00 | 8.01 | 0.07 | 5.77 | 0.04 |
| 2.00 | 1.00 | 0.25 | 7.70 | 0.06 | 5.83 | 0.00 |
| 1.00 | 1.00 | 0.50 | 8.08 | 0.07 | 5.98 | 0.01 |
| 1.00 | 1.50 | 0.25 | 8.03 | 0.07 | 6.11 | 0.00 |
| 1.00 | 4.00 | 4.00 | 8.99 | 0.09 | 6.29 | 0.42 |
| 4.00 | 1.00 | 4.00 | 8.94 | 0.09 | 6.30 | 0.02 |
| 4.00 | 4.00 | 0.12 | 9.02 | 0.09 | 6.42 | 0.00 |
| 4.00 | 0.50 | 0.50 | 9.07 | 0.08 | 6.46 | 0.00 |
| 0.50 | 4.00 | 0.50 | 8.95 | 0.08 | 6.46 | 0.13 |
| 1.00 | 1.00 | 1.00 | 8.94 | 0.08 | 6.53 | 0.07 |
| 1.00 | 1.50 | 2.00 | 9.30 | 0.09 | 6.70 | 0.18 |
| 1.00 | 1.00 | 0.25 | 9.11 | 0.08 | 6.80 | 0.00 |
| 4.00 | 2.00 | 0.12 | 9.57 | 0.09 | 6.81 | 0.00 |
| 2.00 | 4.00 | 0.12 | 9.51 | 0.09 | 6.83 | 0.00 |
| 4.00 | 1.50 | 0.12 | 9.95 | 0.09 | 7.23 | 0.00 |
| 0.50 | 4.00 | 0.25 | 9.75 | 0.09 | 7.28 | 0.00 |
| 0.50 | 4.00 | 1.00 | 10.51 | 0.10 | 7.28 | 0.32 |
| 2.00 | 2.00 | 0.12 | 10.17 | 0.09 | 7.32 | 0.00 |
| 4.00 | 0.50 | 0.25 | 9.84 | 0.09 | 7.39 | 0.00 |
| 2.00 | 0.50 | 0.50 | 10.27 | 0.09 | 7.44 | 0.00 |
| 0.50 | 2.00 | 0.50 | 10.21 | 0.09 | 7.44 | 0.12 |
| 2.00 | 1.50 | 0.12 | 10.36 | 0.09 | 7.51 | 0.00 |
| 4.00 | 0.50 | 1.00 | 10.83 | 0.11 | 7.59 | 0.00 |
| 1.00 | 4.00 | 0.12 | 10.79 | 0.10 | 7.76 | 0.00 |
| 4.00 | 1.00 | 0.12 | 10.66 | 0.10 | 7.79 | 0.00 |
| 0.50 | 1.50 | 0.50 | 10.66 | 0.10 | 7.80 | 0.11 |
| 2.00 | 0.50 | 0.25 | 10.45 | 0.09 | 7.94 | 0.00 |
| 0.50 | 2.00 | 0.25 | 10.54 | 0.09 | 7.97 | 0.01 |
| 1.00 | 2.00 | 0.12 | 11.32 | 0.10 | 8.24 | 0.00 |
| 2.00 | 1.00 | 0.12 | 11.18 | 0.10 | 8.24 | 0.00 |
| 1.00 | 2.00 | 4.00 | 11.98 | 0.12 | 8.41 | 0.31 |
| 0.50 | 1.50 | 0.25 | 11.16 | 0.09 | 8.43 | 0.01 |
| 1.00 | 1.50 | 0.12 | 11.62 | 0.10 | 8.59 | 0.00 |
| 2.00 | 0.50 | 1.00 | 12.25 | 0.11 | 8.69 | 0.01 |
| 2.00 | 1.00 | 4.00 | 12.20 | 0.12 | 8.71 | 0.07 |
| 0.50 | 2.00 | 1.00 | 12.31 | 0.12 | 8.79 | 0.28 |
| 1.00 | 1.00 | 2.00 | 12.53 | 0.12 | 9.03 | 0.14 |
| 1.00 | 0.50 | 0.25 | 12.15 | 0.10 | 9.18 | 0.00 |
| 1.00 | 0.50 | 0.50 | 12.62 | 0.11 | 9.19 | 0.01 |
| 0.50 | 1.00 | 0.25 | 12.22 | 0.10 | 9.20 | 0.01 |
| 0.50 | 1.00 | 0.50 | 12.50 | 0.11 | 9.31 | 0.10 |
| 1.00 | 1.00 | 0.12 | 12.59 | 0.11 | 9.36 | 0.00 |
| 4.00 | 0.50 | 0.12 | 13.13 | 0.11 | 9.79 | 0.00 |
| 0.50 | 4.00 | 0.12 | 13.18 | 0.11 | 9.91 | 0.00 |
| 0.50 | 1.50 | 1.00 | 13.72 | 0.13 | 9.93 | 0.25 |
| 2.00 | 0.50 | 0.12 | 13.71 | 0.12 | 10.19 | 0.00 |
| 1.00 | 1.50 | 4.00 | 14.61 | 0.14 | 10.33 | 0.26 |
| 0.50 | 2.00 | 0.12 | 13.88 | 0.12 | 10.36 | 0.00 |
| 4.00 | 0.50 | 2.00 | 15.33 | 0.15 | 10.72 | 0.00 |
| 0.50 | 4.00 | 2.00 | 15.43 | 0.15 | 10.72 | 0.50 |
| 0.50 | 1.50 | 0.12 | 14.20 | 0.12 | 10.78 | 0.00 |
| 1.00 | 0.50 | 1.00 | 15.73 | 0.15 | 11.33 | 0.04 |
| 1.00 | 0.50 | 0.12 | 15.08 | 0.13 | 11.34 | 0.00 |
| 0.50 | 1.00 | 1.00 | 15.92 | 0.15 | 11.39 | 0.21 |
| 0.50 | 1.00 | 0.12 | 15.16 | 0.13 | 11.58 | 0.00 |
| 4.00 | 0.25 | 0.25 | 16.81 | 0.16 | 11.86 | 0.00 |
| 4.00 | 4.00 | 0.06 | 16.89 | 0.16 | 11.89 | 0.00 |
| 0.50 | 0.50 | 0.25 | 15.90 | 0.13 | 11.90 | 0.01 |
| 0.25 | 4.00 | 0.25 | 17.15 | 0.16 | 12.24 | 0.14 |
| 2.00 | 4.00 | 0.06 | 17.58 | 0.17 | 12.42 | 0.00 |
| 4.00 | 2.00 | 0.06 | 17.41 | 0.17 | 12.46 | 0.00 |
| 2.00 | 2.00 | 0.06 | 18.17 | 0.17 | 12.58 | 0.00 |
| 2.00 | 0.50 | 2.00 | 17.94 | 0.17 | 12.74 | 0.02 |
| 0.50 | 2.00 | 2.00 | 18.08 | 0.18 | 12.83 | 0.41 |
| 0.25 | 2.00 | 0.25 | 17.82 | 0.16 | 12.91 | 0.13 |
| 4.00 | 1.50 | 0.06 | 17.91 | 0.17 | 12.92 | 0.00 |
| 0.50 | 0.50 | 0.50 | 17.91 | 0.16 | 13.02 | 0.07 |
| 2.00 | 1.50 | 0.06 | 18.21 | 0.17 | 13.12 | 0.00 |
| 2.00 | 0.25 | 0.25 | 18.05 | 0.17 | 13.21 | 0.00 |
| 1.00 | 4.00 | 0.06 | 18.72 | 0.17 | 13.28 | 0.00 |
| 4.00 | 1.00 | 0.06 | 18.38 | 0.17 | 13.48 | 0.00 |
| 4.00 | 0.25 | 0.50 | 19.67 | 0.20 | 13.49 | 0.00 |
| 0.25 | 4.00 | 0.50 | 19.33 | 0.19 | 13.50 | 0.35 |
| 0.50 | 0.50 | 0.12 | 18.08 | 0.15 | 13.60 | 0.00 |
| 0.25 | 1.50 | 0.25 | 18.68 | 0.17 | 13.64 | 0.13 |
| 2.00 | 1.00 | 0.06 | 19.19 | 0.18 | 13.71 | 0.00 |
| 1.00 | 2.00 | 0.06 | 19.03 | 0.18 | 13.72 | 0.00 |
| 4.00 | 0.25 | 0.12 | 18.63 | 0.16 | 13.92 | 0.00 |
| 0.50 | 1.50 | 2.00 | 20.15 | 0.20 | 13.99 | 0.37 |
| 1.00 | 1.50 | 0.06 | 19.64 | 0.18 | 14.05 | 0.00 |
| 0.25 | 4.00 | 0.12 | 18.61 | 0.16 | 14.08 | 0.01 |
| 1.00 | 1.00 | 4.00 | 20.27 | 0.19 | 14.18 | 0.19 |
| 1.00 | 1.00 | 0.06 | 19.82 | 0.18 | 14.40 | 0.00 |
| 0.25 | 2.00 | 0.12 | 19.52 | 0.17 | 14.52 | 0.01 |
| 2.00 | 0.25 | 0.12 | 19.57 | 0.17 | 14.54 | 0.00 |
| 1.00 | 0.25 | 0.25 | 20.17 | 0.18 | 14.56 | 0.00 |
| 2.00 | 0.25 | 0.50 | 20.80 | 0.20 | 14.58 | 0.00 |
| 0.25 | 2.00 | 0.50 | 21.01 | 0.20 | 14.80 | 0.32 |
| 0.25 | 1.00 | 0.25 | 20.15 | 0.18 | 14.80 | 0.12 |
| 0.50 | 4.00 | 0.06 | 20.79 | 0.19 | 15.13 | 0.00 |
| 0.25 | 1.50 | 0.12 | 20.06 | 0.17 | 15.21 | 0.01 |
| 4.00 | 0.50 | 0.06 | 20.97 | 0.19 | 15.29 | 0.00 |
| 0.50 | 2.00 | 0.06 | 21.18 | 0.19 | 15.59 | 0.00 |
| 0.25 | 1.50 | 0.50 | 22.21 | 0.21 | 15.63 | 0.30 |
| 2.00 | 0.50 | 0.06 | 21.46 | 0.19 | 15.76 | 0.00 |
| 0.25 | 1.00 | 0.12 | 21.09 | 0.18 | 15.86 | 0.01 |
| 1.00 | 0.25 | 0.12 | 21.34 | 0.18 | 16.01 | 0.00 |
| 0.50 | 1.50 | 0.06 | 21.97 | 0.19 | 16.06 | 0.00 |
| 0.50 | 1.00 | 0.06 | 22.59 | 0.20 | 16.34 | 0.00 |
| 1.00 | 0.50 | 0.06 | 23.09 | 0.21 | 16.85 | 0.00 |
| 4.00 | 0.50 | 4.00 | 24.61 | 0.25 | 16.98 | 0.01 |
| 0.50 | 1.00 | 2.00 | 24.00 | 0.24 | 16.98 | 0.30 |
| 0.50 | 4.00 | 4.00 | 25.20 | 0.25 | 17.19 | 0.63 |
| 1.00 | 0.50 | 2.00 | 24.30 | 0.23 | 17.28 | 0.07 |
| 0.25 | 1.00 | 0.50 | 24.64 | 0.24 | 17.51 | 0.28 |
| 1.00 | 0.25 | 0.50 | 24.78 | 0.24 | 17.73 | 0.01 |
| 0.50 | 0.50 | 1.00 | 24.95 | 0.23 | 17.82 | 0.14 |
| 0.50 | 0.25 | 0.25 | 24.95 | 0.22 | 18.30 | 0.01 |
| 0.50 | 0.50 | 0.06 | 25.07 | 0.22 | 18.49 | 0.00 |
| 0.25 | 0.50 | 0.12 | 24.76 | 0.21 | 18.60 | 0.01 |
| 0.50 | 0.25 | 0.12 | 24.81 | 0.21 | 18.69 | 0.00 |
| 0.25 | 4.00 | 0.06 | 25.44 | 0.23 | 18.74 | 0.00 |
| 0.25 | 0.50 | 0.25 | 25.24 | 0.22 | 18.83 | 0.10 |
| 4.00 | 0.25 | 0.06 | 25.20 | 0.22 | 18.87 | 0.00 |
| 0.25 | 4.00 | 1.00 | 27.71 | 0.28 | 18.98 | 0.55 |
| 4.00 | 0.25 | 1.00 | 27.86 | 0.28 | 19.01 | 0.00 |
| 2.00 | 0.25 | 0.06 | 26.14 | 0.23 | 19.37 | 0.00 |
| 0.25 | 2.00 | 0.06 | 26.01 | 0.22 | 19.39 | 0.00 |
| 0.25 | 1.50 | 0.06 | 26.74 | 0.24 | 19.86 | 0.00 |
| 0.25 | 1.00 | 0.06 | 27.76 | 0.24 | 20.46 | 0.00 |
| 1.00 | 0.25 | 0.06 | 27.31 | 0.24 | 20.53 | 0.00 |
| 0.25 | 2.00 | 1.00 | 29.88 | 0.30 | 20.53 | 0.50 |
| 0.50 | 2.00 | 4.00 | 30.29 | 0.30 | 20.84 | 0.52 |
| 2.00 | 0.25 | 1.00 | 30.26 | 0.30 | 21.34 | 0.00 |
| 2.00 | 0.50 | 4.00 | 30.57 | 0.30 | 21.52 | 0.03 |
| 0.50 | 0.25 | 0.50 | 31.83 | 0.30 | 22.40 | 0.04 |
| 0.25 | 0.50 | 0.06 | 30.35 | 0.26 | 22.63 | 0.00 |
| 0.25 | 1.50 | 1.00 | 32.43 | 0.32 | 22.67 | 0.47 |
| 0.50 | 0.25 | 0.06 | 30.30 | 0.26 | 22.75 | 0.00 |
| 0.25 | 0.50 | 0.50 | 32.20 | 0.31 | 22.92 | 0.22 |
| 4.00 | 0.12 | 0.12 | 32.68 | 0.30 | 23.21 | 0.00 |
| 2.00 | 0.12 | 0.12 | 33.44 | 0.31 | 23.63 | 0.00 |
| 0.25 | 0.25 | 0.12 | 31.53 | 0.27 | 23.65 | 0.01 |
| 0.12 | 4.00 | 0.12 | 33.11 | 0.31 | 23.87 | 0.14 |
| 0.12 | 2.00 | 0.12 | 34.07 | 0.32 | 24.26 | 0.14 |
| 0.50 | 1.50 | 4.00 | 34.88 | 0.34 | 24.48 | 0.46 |
| 0.12 | 1.50 | 0.12 | 34.58 | 0.31 | 24.81 | 0.13 |
| 0.25 | 1.00 | 1.00 | 35.87 | 0.35 | 25.27 | 0.42 |
| 0.12 | 4.00 | 0.25 | 37.00 | 0.37 | 25.61 | 0.36 |
| 1.00 | 0.25 | 1.00 | 36.02 | 0.35 | 25.64 | 0.02 |
| 4.00 | 0.12 | 0.25 | 37.14 | 0.37 | 25.69 | 0.00 |
| 0.25 | 0.25 | 0.25 | 35.72 | 0.32 | 26.14 | 0.07 |
| 0.12 | 1.00 | 0.12 | 36.23 | 0.33 | 26.18 | 0.12 |
| 1.00 | 0.12 | 0.12 | 36.27 | 0.33 | 26.34 | 0.00 |
| 2.00 | 0.12 | 0.25 | 38.84 | 0.39 | 27.08 | 0.00 |
| 0.12 | 4.00 | 0.06 | 36.45 | 0.32 | 27.20 | 0.00 |
| 0.25 | 0.25 | 0.06 | 36.46 | 0.31 | 27.35 | 0.00 |
| 0.12 | 1.50 | 0.25 | 39.11 | 0.38 | 27.59 | 0.33 |
| 4.00 | 0.12 | 0.06 | 36.93 | 0.32 | 27.60 | 0.00 |
| 0.12 | 2.00 | 0.25 | 39.79 | 0.39 | 27.68 | 0.35 |
| 2.00 | 0.12 | 0.06 | 37.30 | 0.32 | 27.78 | 0.00 |
| 0.12 | 2.00 | 0.06 | 37.65 | 0.33 | 27.89 | 0.00 |
| 0.12 | 1.50 | 0.06 | 38.00 | 0.34 | 27.97 | 0.01 |
| 0.50 | 0.50 | 2.00 | 40.64 | 0.40 | 28.38 | 0.19 |
| 0.12 | 0.50 | 0.12 | 40.65 | 0.37 | 28.73 | 0.12 |
| 1.00 | 0.12 | 0.06 | 38.85 | 0.33 | 28.89 | 0.00 |
| 0.12 | 1.00 | 0.06 | 39.12 | 0.34 | 29.13 | 0.01 |
| 0.50 | 0.12 | 0.12 | 40.27 | 0.36 | 29.36 | 0.00 |
| 1.00 | 0.50 | 4.00 | 42.06 | 0.41 | 29.62 | 0.09 |
| 0.50 | 1.00 | 4.00 | 42.36 | 0.42 | 29.63 | 0.37 |
| 1.00 | 0.12 | 0.25 | 42.85 | 0.41 | 29.84 | 0.00 |
| 4.00 | 0.25 | 2.00 | 44.50 | 0.45 | 30.61 | 0.00 |
| 0.12 | 1.00 | 0.25 | 42.87 | 0.41 | 30.77 | 0.33 |
| 0.50 | 0.12 | 0.06 | 42.71 | 0.37 | 32.05 | 0.00 |
| 0.25 | 4.00 | 2.00 | 45.99 | 0.46 | 32.12 | 0.70 |
| 0.12 | 0.50 | 0.06 | 42.40 | 0.36 | 32.18 | 0.01 |
| 0.25 | 0.50 | 1.00 | 48.65 | 0.48 | 33.72 | 0.31 |
| 2.00 | 0.25 | 2.00 | 49.84 | 0.50 | 34.71 | 0.01 |
| 0.50 | 0.12 | 0.25 | 48.94 | 0.47 | 34.73 | 0.01 |
| 0.12 | 0.50 | 0.25 | 49.07 | 0.47 | 34.88 | 0.27 |
| 0.50 | 0.25 | 1.00 | 49.03 | 0.47 | 34.91 | 0.07 |
| 0.25 | 2.00 | 2.00 | 49.97 | 0.50 | 35.23 | 0.63 |
| 0.12 | 4.00 | 0.50 | 52.65 | 0.53 | 35.90 | 0.58 |
| 0.12 | 0.25 | 0.12 | 49.69 | 0.44 | 36.04 | 0.10 |
| 0.25 | 0.25 | 0.50 | 50.69 | 0.47 | 36.10 | 0.14 |
| 0.25 | 0.12 | 0.12 | 50.26 | 0.45 | 36.43 | 0.01 |
| 4.00 | 0.12 | 0.50 | 52.82 | 0.53 | 36.54 | 0.00 |
| 0.25 | 0.12 | 0.06 | 49.08 | 0.42 | 36.82 | 0.00 |
| 0.25 | 1.50 | 2.00 | 53.16 | 0.53 | 36.83 | 0.59 |
| 0.12 | 0.25 | 0.06 | 48.32 | 0.40 | 37.19 | 0.01 |
| 2.00 | 0.12 | 0.50 | 54.47 | 0.55 | 37.69 | 0.00 |
| 0.12 | 2.00 | 0.50 | 54.80 | 0.55 | 37.91 | 0.55 |
| 0.12 | 1.50 | 0.50 | 57.44 | 0.57 | 40.61 | 0.53 |
| 1.00 | 0.12 | 0.50 | 60.27 | 0.61 | 40.92 | 0.00 |
| 1.00 | 0.25 | 2.00 | 59.82 | 0.60 | 41.49 | 0.03 |
| 0.12 | 1.00 | 0.50 | 60.70 | 0.60 | 41.99 | 0.50 |
| 0.25 | 1.00 | 2.00 | 60.33 | 0.59 | 42.22 | 0.52 |
| 0.12 | 0.25 | 0.25 | 63.51 | 0.60 | 44.84 | 0.21 |
| 0.25 | 0.12 | 0.25 | 63.70 | 0.58 | 46.65 | 0.04 |
| 0.12 | 0.12 | 0.06 | 63.47 | 0.54 | 47.17 | 0.01 |
| 0.50 | 0.12 | 0.50 | 72.04 | 0.70 | 50.47 | 0.02 |
| 0.50 | 0.50 | 4.00 | 72.69 | 0.72 | 50.66 | 0.22 |
| 0.12 | 0.50 | 0.50 | 72.76 | 0.70 | 50.89 | 0.42 |
| 0.12 | 0.12 | 0.12 | 71.42 | 0.65 | 52.70 | 0.07 |
| 0.25 | 4.00 | 4.00 | 79.83 | 0.80 | 55.53 | 0.78 |
| 0.25 | 0.25 | 1.00 | 79.37 | 0.75 | 57.01 | 0.19 |
| 4.00 | 0.25 | 4.00 | 81.27 | 0.82 | 57.02 | 0.00 |
| 0.50 | 0.25 | 2.00 | 84.39 | 0.84 | 57.62 | 0.09 |
| 4.00 | 0.12 | 1.00 | 85.18 | 0.86 | 58.77 | 0.00 |
| 0.12 | 4.00 | 1.00 | 86.63 | 0.88 | 59.63 | 0.73 |
| 0.25 | 0.50 | 2.00 | 85.69 | 0.84 | 60.37 | 0.37 |
| 0.25 | 2.00 | 4.00 | 88.57 | 0.89 | 61.06 | 0.70 |
| 2.00 | 0.12 | 1.00 | 90.23 | 0.90 | 62.17 | 0.00 |
| 2.00 | 0.25 | 4.00 | 89.05 | 0.89 | 62.19 | 0.01 |
| 0.12 | 2.00 | 1.00 | 91.40 | 0.91 | 63.45 | 0.70 |
| 0.12 | 1.50 | 1.00 | 93.05 | 0.95 | 64.15 | 0.67 |
| 0.25 | 1.50 | 4.00 | 97.06 | 0.96 | 68.00 | 0.66 |
| 0.12 | 1.00 | 1.00 | 101.14 | 1.00 | 69.28 | 0.63 |
| 0.12 | 0.25 | 0.50 | 96.90 | 0.92 | 69.58 | 0.30 |
| 1.00 | 0.12 | 1.00 | 100.05 | 0.99 | 70.18 | 0.01 |
| 0.25 | 0.12 | 0.50 | 99.49 | 0.95 | 70.42 | 0.07 |
| 0.12 | 0.12 | 0.25 | 99.24 | 0.93 | 71.17 | 0.14 |
| 1.00 | 0.25 | 4.00 | 110.45 | 1.09 | 76.12 | 0.04 |
| 0.25 | 1.00 | 4.00 | 111.00 | 1.10 | 77.52 | 0.57 |
| 0.50 | 0.12 | 1.00 | 120.40 | 1.20 | 83.33 | 0.03 |
| 0.12 | 0.50 | 1.00 | 119.58 | 1.17 | 83.71 | 0.52 |
| 0.25 | 0.25 | 2.00 | 143.75 | 1.42 | 99.77 | 0.22 |
| 4.00 | 0.12 | 2.00 | 150.81 | 1.53 | 103.00 | 0.00 |
| 0.12 | 4.00 | 2.00 | 153.99 | 1.59 | 104.75 | 0.83 |
| 0.50 | 0.25 | 4.00 | 155.15 | 1.52 | 107.42 | 0.10 |
| 0.25 | 0.50 | 4.00 | 156.78 | 1.52 | 111.05 | 0.41 |
| 2.00 | 0.12 | 2.00 | 161.99 | 1.63 | 112.87 | 0.00 |
| 0.12 | 2.00 | 2.00 | 162.69 | 1.62 | 113.27 | 0.78 |
| 0.12 | 0.12 | 0.50 | 162.55 | 1.56 | 113.82 | 0.18 |
| 0.12 | 1.50 | 2.00 | 168.73 | 1.70 | 116.88 | 0.76 |
| 0.12 | 0.25 | 1.00 | 168.58 | 1.66 | 117.31 | 0.37 |
| 0.25 | 0.12 | 1.00 | 168.89 | 1.63 | 119.33 | 0.09 |
| 0.12 | 1.00 | 2.00 | 180.69 | 1.80 | 124.41 | 0.70 |
| 1.00 | 0.12 | 2.00 | 181.32 | 1.81 | 125.22 | 0.01 |
| 0.50 | 0.12 | 2.00 | 219.76 | 2.14 | 153.85 | 0.03 |
| 0.12 | 0.50 | 2.00 | 223.53 | 2.22 | 155.90 | 0.58 |
| 0.25 | 0.25 | 4.00 | 268.27 | 2.59 | 186.44 | 0.23 |
| 4.00 | 0.12 | 4.00 | 287.62 | 2.85 | 196.50 | 0.00 |
| 0.12 | 0.12 | 1.00 | 288.50 | 2.88 | 200.44 | 0.22 |
| 0.12 | 4.00 | 4.00 | 285.97 | 2.76 | 200.90 | 0.88 |
| 2.00 | 0.12 | 4.00 | 307.82 | 2.99 | 214.76 | 0.00 |
| 0.25 | 0.12 | 2.00 | 313.38 | 3.09 | 217.68 | 0.10 |
| 0.12 | 0.25 | 2.00 | 310.65 | 3.05 | 218.58 | 0.41 |
| 0.12 | 1.50 | 4.00 | 314.75 | 3.04 | 219.10 | 0.80 |
| 0.12 | 2.00 | 4.00 | 310.70 | 2.97 | 221.04 | 0.83 |
| 1.00 | 0.12 | 4.00 | 335.68 | 3.23 | 232.55 | 0.01 |
| 0.12 | 1.00 | 4.00 | 332.85 | 3.21 | 235.30 | 0.75 |
| 0.50 | 0.12 | 4.00 | 408.62 | 3.86 | 288.75 | 0.04 |
| 0.12 | 0.50 | 4.00 | 413.18 | 3.90 | 289.44 | 0.61 |
| 0.12 | 0.12 | 2.00 | 537.48 | 5.15 | 384.14 | 0.23 |
| 0.25 | 0.12 | 4.00 | 581.83 | 5.26 | 412.59 | 0.11 |
| 0.12 | 0.25 | 4.00 | 588.86 | 5.34 | 421.26 | 0.42 |
| 0.12 | 0.12 | 4.00 | 982.43 | 8.19 | 739.12 | 0.24 |
Let’s plot the dependence of the median time against each parameter. These plots show a projection of all data onto one parameter dimension.
Since we use quite wide range in each parameter, these plots are very scattered (vertically). This is because for one fixed value of a parameter (e.g. detachment rate of 4 in the right-hand panel) there is a large spread of simulation times due to the range of the other two parameters. These plots are not very informative, let’s try something different.
Here is a plot of median time as a function of formation rate and conversion rate, for two selected values of the detachment rate.
Note that conversion and formation are symmetric, the plots look the same if we swap them.
This is because conversion always follows formation and what matters is the total duration of these two. Since our parameters are rates, mean formation and conversion times are \(1/R_\mathrm{form}\) and \(1/R_\mathrm{conv}\), respectively. Below, we plot the sum of both quantities on the x-axis and the median simulation time on the y-axis. Colour correspond to different detachment rates.
For quick error correction we need high formation and conversion rates, and not very low detachment rate.
Here is a heatmap of formation rate versus conversion rate, for fixed detachment rates of 0.125 and 1. The colour scale and the numbers in boxes show the median simulation time.
The next heatmap shows conversion rate fixed at our best estimate of ~2 per min.
Here we plot distribution of individual simulation times for several sets of parameters. The vertical orange line indicates the median of the distribution.
Here we create a plot for the paper. Following our previous simulations and Vanya’s paper we fix \(R_\mathrm{form} = 2~\mathrm{min}^{-1}\) and \(R_\mathrm{conv} = 2~\mathrm{min}^{-1}\). Shuyu’s observations suggest median bi-orientation time of about 1 min. Hence, for these two fixed parameters we seek the detachment rate with the median run-time closest to 1 min.
The best parameter is \(R_\mathrm{det} = 1~\mathrm{min}^{-1}\). Here is the distribution of individual run times for this parameter. The orange solid line shows the median of the distribution, while the blue dotted lines indicate 95% confidence intervals on the median bi-orientation time from Shuyu’s observations.
## TableGrob (1 x 2) "arrange": 2 grobs
## z cells name grob
## 1 1 (1-1,1-1) arrange gtable[layout]
## 2 2 (1-1,2-2) arrange gtable[layout]
Here is a link to the figure: PDF.
Here is the grid of parameters for model 2:
| replacement rate | 0.125, 0.25, 0.5, 1, 2, 4 |
| knockoff rate | 0.5, 1, 2, 10, 1e+16 |
| conversion rate | 0.125, 0.25, 0.5, 1, 1.5, 2, 4 |
There are 210 combinations of the three parameters.
| replacement | knockoff | conversion | mean.time | se.time | median.time |
|---|---|---|---|---|---|
| 4.00 | 1e+16 | 4.00 | 1.02 | 0.01 | 0.73 |
| 2.00 | 1e+16 | 4.00 | 1.12 | 0.01 | 0.83 |
| 4.00 | 1e+01 | 4.00 | 1.45 | 0.01 | 1.06 |
| 2.00 | 1e+01 | 4.00 | 1.45 | 0.01 | 1.09 |
| 1.00 | 1e+16 | 4.00 | 1.57 | 0.01 | 1.18 |
| 1.00 | 1e+01 | 4.00 | 1.81 | 0.02 | 1.34 |
| 2.00 | 1e+16 | 2.00 | 1.99 | 0.02 | 1.41 |
| 4.00 | 1e+16 | 2.00 | 2.25 | 0.02 | 1.53 |
| 1.00 | 1e+16 | 2.00 | 2.24 | 0.02 | 1.67 |
| 2.00 | 1e+01 | 2.00 | 2.41 | 0.02 | 1.75 |
| 0.50 | 1e+16 | 4.00 | 2.51 | 0.02 | 1.80 |
| 1.00 | 1e+01 | 2.00 | 2.55 | 0.02 | 1.92 |
| 2.00 | 1e+16 | 1.50 | 2.76 | 0.03 | 1.93 |
| 4.00 | 1e+01 | 2.00 | 2.88 | 0.03 | 2.04 |
| 0.50 | 1e+01 | 4.00 | 2.79 | 0.02 | 2.08 |
| 1.00 | 1e+16 | 1.50 | 2.81 | 0.02 | 2.09 |
| 2.00 | 1e+01 | 1.50 | 3.14 | 0.03 | 2.22 |
| 4.00 | 1e+16 | 1.50 | 3.35 | 0.03 | 2.27 |
| 1.00 | 1e+01 | 1.50 | 3.10 | 0.03 | 2.28 |
| 0.50 | 1e+16 | 2.00 | 3.09 | 0.03 | 2.30 |
| 1.00 | 2e+00 | 4.00 | 3.04 | 0.02 | 2.35 |
| 2.00 | 2e+00 | 4.00 | 3.15 | 0.03 | 2.35 |
| 0.50 | 1e+01 | 2.00 | 3.37 | 0.03 | 2.52 |
| 0.50 | 1e+16 | 1.50 | 3.53 | 0.03 | 2.68 |
| 0.50 | 2e+00 | 4.00 | 3.75 | 0.03 | 2.80 |
| 4.00 | 2e+00 | 4.00 | 3.99 | 0.04 | 2.82 |
| 0.50 | 1e+01 | 1.50 | 3.84 | 0.03 | 2.86 |
| 1.00 | 1e+16 | 1.00 | 4.04 | 0.04 | 2.87 |
| 4.00 | 1e+01 | 1.50 | 4.20 | 0.04 | 2.96 |
| 1.00 | 2e+00 | 2.00 | 3.98 | 0.03 | 2.97 |
| 2.00 | 1e+16 | 1.00 | 4.45 | 0.04 | 3.06 |
| 1.00 | 1e+01 | 1.00 | 4.35 | 0.04 | 3.15 |
| 0.25 | 1e+16 | 4.00 | 4.58 | 0.04 | 3.27 |
| 2.00 | 2e+00 | 2.00 | 4.51 | 0.04 | 3.30 |
| 0.50 | 1e+16 | 1.00 | 4.51 | 0.04 | 3.34 |
| 0.25 | 1e+01 | 4.00 | 4.71 | 0.04 | 3.39 |
| 0.50 | 2e+00 | 2.00 | 4.52 | 0.04 | 3.42 |
| 1.00 | 2e+00 | 1.50 | 4.73 | 0.04 | 3.53 |
| 2.00 | 1e+01 | 1.00 | 5.07 | 0.05 | 3.55 |
| 0.50 | 1e+01 | 1.00 | 4.81 | 0.04 | 3.58 |
| 0.25 | 1e+16 | 2.00 | 5.08 | 0.05 | 3.67 |
| 1.00 | 1e+00 | 4.00 | 5.11 | 0.05 | 3.71 |
| 0.50 | 2e+00 | 1.50 | 4.99 | 0.04 | 3.78 |
| 0.25 | 1e+01 | 2.00 | 5.30 | 0.05 | 3.82 |
| 0.50 | 1e+00 | 4.00 | 5.28 | 0.05 | 3.98 |
| 0.25 | 1e+16 | 1.50 | 5.47 | 0.05 | 4.03 |
| 2.00 | 2e+00 | 1.50 | 5.67 | 0.05 | 4.13 |
| 0.25 | 1e+01 | 1.50 | 5.67 | 0.05 | 4.15 |
| 0.25 | 2e+00 | 4.00 | 5.71 | 0.05 | 4.22 |
| 4.00 | 1e+16 | 1.00 | 6.31 | 0.06 | 4.26 |
| 2.00 | 1e+00 | 4.00 | 6.15 | 0.06 | 4.37 |
| 4.00 | 2e+00 | 2.00 | 6.16 | 0.06 | 4.41 |
| 1.00 | 2e+00 | 1.00 | 6.14 | 0.05 | 4.54 |
| 0.50 | 1e+00 | 2.00 | 6.11 | 0.05 | 4.59 |
| 1.00 | 1e+00 | 2.00 | 6.19 | 0.06 | 4.59 |
| 0.25 | 2e+00 | 2.00 | 6.17 | 0.05 | 4.61 |
| 0.50 | 2e+00 | 1.00 | 6.15 | 0.05 | 4.67 |
| 0.25 | 1e+16 | 1.00 | 6.38 | 0.06 | 4.77 |
| 0.25 | 1e+01 | 1.00 | 6.48 | 0.06 | 4.82 |
| 0.50 | 1e+00 | 1.50 | 6.67 | 0.06 | 4.96 |
| 0.25 | 2e+00 | 1.50 | 6.72 | 0.06 | 5.00 |
| 4.00 | 1e+01 | 1.00 | 7.38 | 0.07 | 5.06 |
| 1.00 | 1e+00 | 1.50 | 7.11 | 0.06 | 5.24 |
| 0.25 | 1e+00 | 4.00 | 6.98 | 0.06 | 5.24 |
| 2.00 | 2e+00 | 1.00 | 7.96 | 0.08 | 5.58 |
| 0.25 | 2e+00 | 1.00 | 7.55 | 0.06 | 5.66 |
| 4.00 | 2e+00 | 1.50 | 8.05 | 0.08 | 5.71 |
| 0.25 | 1e+00 | 2.00 | 7.54 | 0.06 | 5.73 |
| 0.50 | 1e+16 | 0.50 | 7.94 | 0.07 | 5.76 |
| 2.00 | 1e+00 | 2.00 | 8.17 | 0.08 | 5.89 |
| 0.12 | 1e+16 | 4.00 | 8.50 | 0.08 | 6.00 |
| 0.50 | 1e+00 | 1.00 | 7.99 | 0.07 | 6.01 |
| 0.12 | 1e+01 | 4.00 | 8.62 | 0.08 | 6.05 |
| 0.25 | 1e+00 | 1.50 | 8.11 | 0.07 | 6.07 |
| 0.50 | 1e+01 | 0.50 | 8.43 | 0.08 | 6.14 |
| 1.00 | 1e+16 | 0.50 | 9.16 | 0.09 | 6.29 |
| 4.00 | 1e+00 | 4.00 | 8.97 | 0.09 | 6.29 |
| 0.12 | 1e+16 | 2.00 | 9.04 | 0.09 | 6.42 |
| 0.12 | 1e+01 | 2.00 | 9.17 | 0.09 | 6.50 |
| 0.50 | 5e-01 | 4.00 | 8.98 | 0.08 | 6.52 |
| 1.00 | 1e+00 | 1.00 | 8.96 | 0.08 | 6.57 |
| 0.25 | 1e+16 | 0.50 | 8.93 | 0.08 | 6.69 |
| 2.00 | 1e+00 | 1.50 | 9.46 | 0.09 | 6.72 |
| 0.25 | 1e+00 | 1.00 | 8.99 | 0.08 | 6.73 |
| 0.12 | 1e+16 | 1.50 | 9.37 | 0.09 | 6.74 |
| 1.00 | 1e+01 | 0.50 | 9.58 | 0.09 | 6.83 |
| 0.25 | 1e+01 | 0.50 | 9.35 | 0.08 | 6.91 |
| 0.12 | 1e+01 | 1.50 | 9.63 | 0.09 | 6.92 |
| 0.12 | 2e+00 | 4.00 | 9.57 | 0.09 | 7.00 |
| 0.12 | 1e+16 | 1.00 | 9.93 | 0.09 | 7.18 |
| 0.25 | 5e-01 | 4.00 | 9.73 | 0.08 | 7.21 |
| 0.50 | 5e-01 | 2.00 | 10.19 | 0.09 | 7.41 |
| 0.12 | 2e+00 | 2.00 | 10.24 | 0.09 | 7.43 |
| 1.00 | 5e-01 | 4.00 | 10.48 | 0.10 | 7.45 |
| 0.50 | 2e+00 | 0.50 | 10.16 | 0.09 | 7.49 |
| 0.12 | 1e+01 | 1.00 | 10.35 | 0.09 | 7.58 |
| 0.12 | 2e+00 | 1.50 | 10.45 | 0.09 | 7.67 |
| 0.12 | 1e+00 | 4.00 | 10.61 | 0.09 | 7.80 |
| 0.50 | 5e-01 | 1.50 | 10.85 | 0.10 | 7.84 |
| 0.25 | 2e+00 | 0.50 | 10.49 | 0.09 | 7.84 |
| 0.25 | 5e-01 | 2.00 | 10.69 | 0.09 | 8.06 |
| 0.12 | 2e+00 | 1.00 | 11.13 | 0.10 | 8.20 |
| 0.25 | 5e-01 | 1.50 | 11.08 | 0.10 | 8.28 |
| 0.12 | 1e+00 | 2.00 | 11.31 | 0.10 | 8.30 |
| 4.00 | 1e+00 | 2.00 | 12.04 | 0.12 | 8.43 |
| 0.12 | 1e+00 | 1.50 | 11.54 | 0.10 | 8.48 |
| 2.00 | 1e+16 | 0.50 | 12.49 | 0.13 | 8.54 |
| 4.00 | 2e+00 | 1.00 | 12.19 | 0.12 | 8.61 |
| 1.00 | 2e+00 | 0.50 | 12.17 | 0.12 | 8.68 |
| 1.00 | 5e-01 | 2.00 | 12.09 | 0.11 | 8.74 |
| 2.00 | 1e+00 | 1.00 | 12.46 | 0.12 | 8.89 |
| 0.50 | 1e+00 | 0.50 | 12.29 | 0.11 | 8.99 |
| 0.12 | 1e+00 | 1.00 | 12.28 | 0.11 | 9.12 |
| 0.25 | 5e-01 | 1.00 | 12.25 | 0.10 | 9.14 |
| 0.50 | 5e-01 | 1.00 | 12.51 | 0.11 | 9.17 |
| 0.12 | 1e+16 | 0.50 | 12.46 | 0.11 | 9.20 |
| 0.25 | 1e+00 | 0.50 | 12.26 | 0.11 | 9.30 |
| 2.00 | 1e+01 | 0.50 | 13.62 | 0.14 | 9.46 |
| 0.12 | 1e+01 | 0.50 | 12.59 | 0.11 | 9.56 |
| 1.00 | 5e-01 | 1.50 | 13.50 | 0.13 | 9.70 |
| 0.12 | 5e-01 | 4.00 | 13.22 | 0.11 | 9.88 |
| 4.00 | 1e+00 | 1.50 | 14.58 | 0.14 | 10.21 |
| 0.12 | 2e+00 | 0.50 | 13.71 | 0.12 | 10.22 |
| 0.12 | 5e-01 | 2.00 | 13.74 | 0.12 | 10.28 |
| 2.00 | 5e-01 | 4.00 | 15.10 | 0.15 | 10.62 |
| 0.12 | 5e-01 | 1.50 | 14.29 | 0.12 | 10.75 |
| 1.00 | 1e+00 | 0.50 | 15.93 | 0.15 | 11.33 |
| 1.00 | 5e-01 | 1.00 | 15.95 | 0.15 | 11.39 |
| 0.12 | 1e+00 | 0.50 | 15.17 | 0.13 | 11.47 |
| 0.25 | 1e+16 | 0.25 | 16.12 | 0.15 | 11.48 |
| 0.12 | 5e-01 | 1.00 | 15.29 | 0.13 | 11.56 |
| 0.25 | 1e+01 | 0.25 | 16.54 | 0.15 | 11.88 |
| 0.25 | 5e-01 | 0.50 | 16.06 | 0.14 | 11.92 |
| 0.50 | 1e+16 | 0.25 | 17.61 | 0.18 | 12.17 |
| 2.00 | 5e-01 | 2.00 | 18.18 | 0.18 | 12.51 |
| 2.00 | 2e+00 | 0.50 | 18.29 | 0.18 | 12.70 |
| 0.50 | 1e+01 | 0.25 | 18.59 | 0.18 | 12.84 |
| 0.50 | 5e-01 | 0.50 | 17.84 | 0.16 | 12.99 |
| 0.25 | 2e+00 | 0.25 | 17.99 | 0.17 | 13.02 |
| 0.12 | 1e+16 | 0.25 | 17.98 | 0.16 | 13.38 |
| 0.12 | 1e+01 | 0.25 | 18.15 | 0.16 | 13.49 |
| 4.00 | 1e+16 | 0.50 | 20.05 | 0.21 | 13.64 |
| 0.12 | 5e-01 | 0.50 | 18.08 | 0.15 | 13.65 |
| 2.00 | 5e-01 | 1.50 | 19.99 | 0.19 | 14.00 |
| 4.00 | 1e+00 | 1.00 | 20.05 | 0.19 | 14.18 |
| 0.25 | 1e+00 | 0.25 | 20.32 | 0.19 | 14.48 |
| 0.12 | 2e+00 | 0.25 | 19.33 | 0.17 | 14.55 |
| 4.00 | 1e+01 | 0.50 | 21.48 | 0.22 | 14.89 |
| 0.50 | 2e+00 | 0.25 | 21.38 | 0.21 | 15.00 |
| 0.12 | 1e+00 | 0.25 | 21.14 | 0.18 | 16.00 |
| 1.00 | 1e+16 | 0.25 | 24.71 | 0.25 | 16.79 |
| 2.00 | 5e-01 | 1.00 | 24.15 | 0.23 | 17.13 |
| 2.00 | 1e+00 | 0.50 | 24.70 | 0.24 | 17.29 |
| 0.50 | 1e+00 | 0.25 | 24.48 | 0.23 | 17.34 |
| 1.00 | 1e+01 | 0.25 | 25.42 | 0.26 | 17.45 |
| 4.00 | 5e-01 | 4.00 | 24.79 | 0.25 | 17.49 |
| 0.12 | 5e-01 | 0.25 | 24.18 | 0.20 | 18.12 |
| 1.00 | 5e-01 | 0.50 | 25.26 | 0.24 | 18.15 |
| 0.25 | 5e-01 | 0.25 | 25.19 | 0.22 | 18.51 |
| 4.00 | 2e+00 | 0.50 | 30.06 | 0.30 | 20.90 |
| 1.00 | 2e+00 | 0.25 | 30.15 | 0.30 | 21.06 |
| 4.00 | 5e-01 | 2.00 | 30.78 | 0.30 | 21.81 |
| 0.12 | 1e+16 | 0.12 | 31.84 | 0.30 | 22.72 |
| 0.50 | 5e-01 | 0.25 | 31.75 | 0.29 | 23.02 |
| 0.12 | 1e+01 | 0.12 | 32.57 | 0.31 | 23.24 |
| 4.00 | 5e-01 | 1.50 | 33.96 | 0.34 | 23.73 |
| 0.12 | 2e+00 | 0.12 | 34.64 | 0.32 | 24.77 |
| 0.25 | 1e+16 | 0.12 | 36.08 | 0.37 | 24.79 |
| 0.25 | 1e+01 | 0.12 | 36.05 | 0.36 | 24.81 |
| 1.00 | 1e+00 | 0.25 | 35.62 | 0.35 | 24.97 |
| 0.12 | 1e+00 | 0.12 | 36.07 | 0.33 | 26.15 |
| 0.25 | 2e+00 | 0.12 | 38.46 | 0.38 | 26.45 |
| 2.00 | 1e+16 | 0.25 | 41.07 | 0.43 | 27.33 |
| 2.00 | 1e+01 | 0.25 | 42.35 | 0.43 | 29.18 |
| 2.00 | 5e-01 | 0.50 | 40.67 | 0.39 | 29.33 |
| 4.00 | 1e+00 | 0.50 | 42.60 | 0.43 | 29.44 |
| 0.12 | 5e-01 | 0.12 | 40.94 | 0.37 | 29.49 |
| 4.00 | 5e-01 | 1.00 | 42.44 | 0.42 | 29.84 |
| 0.25 | 1e+00 | 0.12 | 42.63 | 0.41 | 29.93 |
| 0.50 | 1e+16 | 0.12 | 49.46 | 0.51 | 33.41 |
| 2.00 | 2e+00 | 0.25 | 50.56 | 0.51 | 33.74 |
| 1.00 | 5e-01 | 0.25 | 48.47 | 0.46 | 34.09 |
| 0.25 | 5e-01 | 0.12 | 49.05 | 0.47 | 34.66 |
| 0.50 | 1e+01 | 0.12 | 51.68 | 0.53 | 35.65 |
| 0.50 | 2e+00 | 0.12 | 54.52 | 0.55 | 37.39 |
| 2.00 | 1e+00 | 0.25 | 60.48 | 0.60 | 41.48 |
| 0.50 | 1e+00 | 0.12 | 61.45 | 0.61 | 42.82 |
| 4.00 | 5e-01 | 0.50 | 71.02 | 0.71 | 48.93 |
| 4.00 | 1e+16 | 0.25 | 72.14 | 0.74 | 49.76 |
| 4.00 | 1e+01 | 0.25 | 75.97 | 0.77 | 51.51 |
| 0.50 | 5e-01 | 0.12 | 72.91 | 0.72 | 51.74 |
| 1.00 | 1e+16 | 0.12 | 80.87 | 0.83 | 55.29 |
| 1.00 | 1e+01 | 0.12 | 83.77 | 0.85 | 58.23 |
| 2.00 | 5e-01 | 0.25 | 83.71 | 0.81 | 59.25 |
| 1.00 | 2e+00 | 0.12 | 90.68 | 0.91 | 62.90 |
| 4.00 | 2e+00 | 0.25 | 90.33 | 0.89 | 63.25 |
| 1.00 | 1e+00 | 0.12 | 97.74 | 0.98 | 67.48 |
| 4.00 | 1e+00 | 0.25 | 110.52 | 1.10 | 77.20 |
| 1.00 | 5e-01 | 0.12 | 118.63 | 1.17 | 81.62 |
| 2.00 | 1e+16 | 0.12 | 143.87 | 1.47 | 97.36 |
| 2.00 | 1e+01 | 0.12 | 149.35 | 1.54 | 101.89 |
| 4.00 | 5e-01 | 0.25 | 155.35 | 1.51 | 109.51 |
| 2.00 | 2e+00 | 0.12 | 164.40 | 1.66 | 113.47 |
| 2.00 | 1e+00 | 0.12 | 181.06 | 1.81 | 125.70 |
| 2.00 | 5e-01 | 0.12 | 217.40 | 2.12 | 154.23 |
| 4.00 | 1e+16 | 0.12 | 270.70 | 2.68 | 186.34 |
| 4.00 | 1e+01 | 0.12 | 279.42 | 2.74 | 195.13 |
| 4.00 | 2e+00 | 0.12 | 302.45 | 2.91 | 212.57 |
| 4.00 | 1e+00 | 0.12 | 335.98 | 3.25 | 231.54 |
| 4.00 | 5e-01 | 0.12 | 408.97 | 3.85 | 290.47 |
Here is a plot of median time as a function of replacement rate and conversion rate, for two selected values of the knockoff rate.
Here is the heatmap for knockoff rate fixed at 1e16 (no delay) and 1.
Very quick error correction times are achieved for high conversion rates. The knock-off delay has detrimental effect on these times.
The next heatmap shows conversion rate fixed at our best estimate of ~2 per min.
Here we create a plot for the paper. Following our previous simulations and Vanya’s paper we fix \(R_\mathrm{repl} = 2~\mathrm{min}^{-1}\) and \(R_\mathrm{conv} = 2~\mathrm{min}^{-1}\). Shuyu’s observations suggest median bi-orientation time of about 1 min. Hence, for these two fixed parameters we seek the detachment rate with the median run-time closest to 1 min.
The best parameter is \(R_\mathrm{ko} = 10^{16}~\mathrm{min}^{-1}\). Here is the distribution of individual run times for this parameter. The orange solid line shows the median of the distribution, while the blue dotted lines indicate 95% confidence intervals on the median bi-orientation time from Shuyu’s observations.
## TableGrob (1 x 2) "arrange": 2 grobs
## z cells name grob
## 1 1 (1-1,1-1) arrange gtable[layout]
## 2 2 (1-1,2-2) arrange gtable[layout]
Here is a link to the figure: PDF.
In this Section we investigate how long the system spends in configuration where both kinetochores are detached at the same time (model M1 only).
First, we calculate the fraction of the total time spend as fully detached. Here are plots showing the detached fraction across the range of parameters.
Detached periods are getting long for small formation rate and high detachment rate (as expected).
Now, we turn to individual fully detached periods. We notice that the individual detached times depend almost exclusively on the formation rate. We can clearly see this when plotting median detached time:
Hence, we can simplify distribution plots. We fix \(R_\mathrm{conv} = 2\) and \(R_\mathrm{det} = 1~\mathrm{min}^{-1}\) and show the full range of formation rate (the same distribution in linear of logarithmic scale).
Next, we find the diffusion range, defined as \(\sigma = \sqrt{6 D t}\), where \(D = 0.1~\mu \mathrm{m}^2\mathrm{min}^{-1}\) is the diffusion coefficient and \(t\) is the detachment time. We fix \(R_\mathrm{conv} = 2\) and \(R_\mathrm{det} = 1~\mathrm{min}^{-1}\) and show the full range of formation rate (the same distribution in linear of logarithmic scale).
Let us consider \(R_\mathrm{form} = 1~\mathrm{min}^{-1}\). For this, the median diffusion range is \(\sigma\) = 0.45 \(\mu\)m. This is the standard deviation of the Gaussian distribution. Hence, about 95% of KTs will not strand farther than about \(2\sigma\) = 0.91 \(\mu\)m.
Here we take the \(\sigma\) distribution from the above figure (for each formation rate). Then, for each individual \(\sigma\) from the distribution we build the corresponding Gaussian distribution, describing the local diffusion distance. We add up all the distributions to create a “pile-up”. This describes the average distribution of diffusion distances across all observed detachment times. Here is the result. The vertical red lines show the size of the nucleus. Again, we fix \(R_\mathrm{conv} = 2\) and \(R_\mathrm{det} = 1~\mathrm{min}^{-1}\)
The table below shows the fraction of KTs that diffuse within the nucleus (that is, the area of each distribution within the red lines):
| Form rate | Fraction in nuc |
|---|---|
| 4 | 1.00 |
| 2 | 0.99 |
| 1 | 0.96 |
| 0.5 | 0.90 |
| 0.25 | 0.80 |
| 0.125 | 0.69 |
The model described in this document simulates error correction for one chromosome (sister chromatids). In yeast, we have 16 chromosomes. What is the total time required to correct all chromosomes and establish proper bi-orientation? If error correction for all of them is done in parallel, then the total time is the longest of the individual times. We can model this by randomly sampling 16 timelines from our simulation results and finding the maximum time. NOTE: this is a very conservative approach assuming that all chromosomes start in an erroneous configuration.
We select data with the following (default) parameters (units are 1 / minute).
These plots show the distribution of the maximum time across 16 randomly selected samples: